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We study spin ^ hadronic states in quenched lattice QCD to search for a possible S' = +1 
pentaquark resonance. Simulations are carried out on 8® x 24, 10® x 24, 12® x 24 and 16® x 24 lattices 
at 13=5.7 at the quenched level with the standard plaquette gauge action and the Wilson quark 
action. We adopt a Dirichlet boundary condition in the time direction for the quark to circumvent 
the possible contaminations due to the (anti)periodic boundary condition for the quark field, which 
are peculiar to the pentaquark. By diagonalizing the 2x2 correlation matrices constructed from 
two independent operators with the quantum numbers (/, J) = (0, |), we successfully obtain the 
energies of the lowest state and the 2nd-lowest state in this channel. The analysis of the volume 
dependence of the energies and spectral weight factors indicates that a resonance state is likely to 
exist slightly above the NK threshold in {I, J^) = (0, | ) channel. 


I. INTRODUCTION 

After the discovery P| of 0'''(154O) followed by the other experiments iiiiiHili , identifying th^roperties 
of the particle is one of the central problems in hadron physics. While the isospin of G'*' is likely to be zero jj, the spin 
and the parity and the origin of its tiny width still remain open questions |inllll| . In spite of many theoretical studies 
on 0+ [iO. 113 ) llJj llJ> Ell Ea 1 the nature of this exotic particle, including the very existence of the particle, is still 
controversial. Among theoretical approaches, the lattice QCD calculation is considered as one of the most reliable 
ab initio methods for studying the properties of hadronic states, which has been very successful in reproducing the 
non-exotic hadron mass spectra Up to now, several lattice QCD studies have been reported, which aim to look 
for pentaquarks in various different ways. However, the conclusions are unfortunately contradictory with each other. 
On one hand, the authors in Refs. 0,[3,|23 conclude that the parity of 0+ is likely to be negati ve, while in Ref. 
the state with the similar mass to G'*' in the positive parity channel is reported. In Refs. [HjEII, the absence of 0+ 
is suggested. 

One of the difficulties in the spectroscopy calculation with lattice QCD arises from the fact that the hadron masses 
suffer from systematic errors due to the discretization, the chiral extrapolation, the quenching effect, the finite volume 
effect and the contaminations from higher excited-states. The difficulty specific to the present problem is that the 
signal of 0'*' is embedded in the discrete spectrum of NK scattering states in finite volume. In order to verify the 
existence of a resonance state, one needs to isolate the first few low energy states including the lowest NK scattering 
state, identify a resonance state and study its volume dependence which can distinguish itself from other scattering 
states. Therefore, ideally one should extract multistates from a high statistics unquenched calculation for several 
different physical volumes, where both the continuum and the chiral limits are taken. However, due to enormous 
computational costs, so far there are no lattice QCD study which performs all these steps. 

In the present situation where even the very existence of the resonance state is theoretically in dispute, the primary 
task is to provide evidences which distinguish the candidate resonance state from a scattering state. As long as other 
systematic errors only affect the numerical values of the masses but not the characteristic evidences of the resonance 
state, they may be neglected. Even so, the isolation of the first few low energy states and the study of the volume 
dependence is a minimum requisite. 

Therefore, at this stage as a first step towards a more complete analysis, we propose to focus only on analyses using 
rather heavy quarks on coarse quenched lattices but with a good statistics. By such a strategy, we can afford taking 
several different lattice volumes with thousands of gauge configurations so that the careful separation of states and 
the studies of volume dependence are possible. Although giving well controlled continuum and chiral extrapolations 
may be important, we simply assume that the contents in spectra would not be drastically changed, although there 
are some cases where level crossings of resonance states occur as the quark masses decrease [^. Even with such a 
compromise, we can hopefully learn about the existence and much of the qualitative properties of 0^. 

In this paper, we study (/, J) = (0, i) channel in quenched lattice QCD to search for possible resonance states. We 
adopt two independent operators with 1 = 0 and J = \ and diagonalize the 2x2 correlation matrices by the variational 
method for all the combinations of lattice sizes and quark masses to extract the 2nd-lowest state slightly above the 
NK threshold in this channel. After the careful separation of the states, we investigate the volume dependence of the 
energy as well as the spectral weight of each state so that we can distinguish the resonance state from simple 
scattering states. 
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The paper is organized as follows. We present the formalism used in the analysis in Sec. inland show the simulation 
conditions in Sec. EH The process of the analysis is shown in Sec. CYI Secs. IVIVIIII are devoted to the interpretation of 
the obtained results and the verification of the existence of a resonance state, as well as some checks on the consistency 
and the reliability of the obtained data. In Sec. HU we discuss the operator dependence of the results and compare 
our results with the previous works. We finally summarize the paper in Sec. HI In Appendix, we show the result of 
another trial to estimate the volume dependences of the spectral weights, which requires no multi-exponential ht. 


II. FORMALISM 

As ©■*■ lies above the NK threshold, any hadron correlators, which have the signal, also contain the discrete-level 
NK scattering states in a finite volume lattice. In order to isolate the resonance state from the scattering states, one 
needs to extract at least two states before anything else. 

Since a double-exponential ht of a single correlator becomes numerically ambiguous, we adopt the variational method 
using correlation matrices constructed from independent operators 1^ 1^ 1^ . A set of independent operators 
>{^snk} sinks andjOgr^;.} for sources, is needed to construct correlation matrices (0)), which 

can be decomposed into the sum over the energy eigenstates |*) as 

Cl AT) = {OUmoilm = E E = {Cl^AiT)C.,,)u, (1) 

^ j 

with the general matrices which depend on the operators as 

cLr = {o\oUA,arcji = (j|o/;jo), (2) 

and the diagonal matrix 

A(T),, ^ (3) 

From the product 

C-AT + l)C(T) = C',-iA(-I)C's.c, (4) 

we can extract the energies {Ei} as the logarithm of eigenvalues {e^’} of the matrix C~^{T -|- l)C(T). 

While there are N independent operators for the correlation matrix, the number of the intermediate states |*) 
which effectively contribute to this matrix may differ from N in general. Let us call this number as Neg. If Ngg is 
larger than N, the higher excited-states are non negligible and their contaminations give rise to a T-dependence of 
eigenvalues as jf qj^ other hand Ngg is smaller than N, C becomes non-invertible so that the extracted 

energies become numerically fairly unstable and we cannot extract all the N eigenvalues. In order to have a reliable 
extraction of states, we therefore need to find an appropriate window of T (Tmin < T < Tmax) so that Wff = N. 
(Of course, even in the case when N > N^s, we can extract Nes eigenvalues with the reduced N^ff x N^s correlation 
matrices.) The stability of against T is expected in this T range and we can obtain N eigenenergies {Ei} 

(0 < j < A — 1) by fitting the eigenvalues e^^T) 9,3 Ei = Ei{T) = ln(e^*(^)) in Tmin < T < Tmax- Since finding 
the stability of the energies against T in noisy data may suffer from uncontrollable biases, the result could be quite 
subjective. In order to avoid such biases, one should impose some concrete criteria to judge the stability as will be 
explained in later sections and select only those data which satisfy the criteria. After the separation of the states, we 
need to distinguish a possible resonance state from NK scattering states by the volume dependence of each state. It 
is expected that the energies of resonance states have small volume dependence, while the energies of NK scattering 

states are expected to scale as ^+ |^np -I- -I- |^np according to the relative momentum between 

N and K on a finite periodic lattice, provided that the NK interaction is weak and negligible which is indeed the case 
for the leading order in chiral perturbation theory. 

Although the variational method is powerful for extracting the energy spectrum, one can obtain only part of the 
information on the spectral weights C. In order to extract the spectral weights, we also perform constrained double 
exponential fits using the energies from the variational method as inputs. 


III. LATTICE SET UP 

We carry out simulations on four different sizes of lattices, 8^ x24, lO^x 24, 12^x 24 and 16^x 24 with 2900, 2900, 1950 
and 950 gauge conhgurations using the standard plaquette (Wilson) gauge action at (3 = 5.7 and the Wilson quark 
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action. The hopping parameters for the quarks are Ks)=( 0.1600, 0.1650), (0.1625,0.1650), (0.1650,0.1650), 
(0.1600,0.1600) and (0.1650,0.1600), which correspond to the current quark masses (m„^(i,ms) ^ (240,100), 
(170,100), (100,100), (240,240) and (100,240), respectively in the unit of MeV The lattice spacing a from 
the Sommer scale is set to be 0.17 fm, which implies the physical lattice sizes are 1.4^ x 4.0 fm"'’, 1.7^ x 4.0 fm"^, 
2.0^ X 4.0 fm^ and 2.7^ x 4.0 fm^. 

We adopt the following two operators used in Ref. 0 for the interpolating operators at the sink 

0^(x) = £°‘^^[u]'(x)Cj 5 db(x)]{ue(x)[s^(x)j 5 dc(x)] - (u ^ d)}, (5) 

which is expected to have a larger overlap with G"'' state, and 

e^(x) = £°’^^[uJ(x)Cj 5 db(x)]{uc(x)[s^(x)j 5 de(x)] - (u ^ d)}, (6) 

which we expect to have larger overlaps with NK scattering states. Here, the Dirac fields u(x), d(x) and s(a;) are up, 
down and strange quark fields, respectively and the Roman alphabets {a,b,c,e} denote color indices. For measuring 
the energy spectrum, the two operators at the source {Ofrl} 9'''® chosen to be and 0waii(O defined using 

spatially spread quark fields g(x) with the Coulomb gauge: 

©wallW = X! s‘'^‘'[uI(xi)Cj5db(x2)]{Ue(X3)[^(x4)75dc(X5)]-(u^d)}, (7) 

and 

0lall(i) = £°‘^‘'{uJ{xi)C75db{x2)]{Uc{X3)[Mx4h5de{x5)] - {u ^ d)}. (8) 

\ / Xi ~X5 


The above operators give a 2 x 2 correlation matrix in the channel with the quantum number of {I,J) = (0, ^). 
We note here that the baryonic correlators have the spinor indices, which we omit in the paper, and they contain 
the propagations of both the positive and negative parity particles. For the parity projection, we simply multiply 
the correlators by i(l ± 70 ) and extract the contributions proportional to (1 ± 70 ) from the negative-parity and 
positive-parity particles, respectively. 

We fix the source operator 0waii(i) on t = tgrc = 6 plane to reduce the effect of the Dirichlet boundary on t = 0 
plane We adopt the operators operators, which is summed over all space to project out 

the zero-momentum states. We finally calculate 

C^^(T) =^(0^(x,T + t,rc)^wall(tsrc)). (9) 

X 

Using two independent operators, we can extract the first two states, namely the lowest and the next-lowest states. 
The lowest state is considered to be the “lowest” NK scattering state. In order to extract a possible resonance state 
with controlled systematic errors, we need to choose the physical volume of the lattice in an appropriate range. If we 
choose L to be too large, the resonance state becomes heavier than the 2nd-lowest NK scattering state whose energy 

is naively expected to scale as ^Mfj + + \JMj. + according to the spatial lattice extent L. In this case 

we need to extract the 3rd state using a 3 x 3 correlation matrix, which requires more computational time. The energy 
difference between the lowest and the next-lowest NK scattering states + \J— {Mn + Mk), 

for example, ranges from 180 MeV to 860 MeV in 1.4fm < L < 3.5fm. Taking into account that B'*' lies about 
100 MeV above the NK threshold, we take 3.5 fm as the upper limit of L. On the other hand, if we choose L too 
small, unwanted finite-volume artifacts from the finite sizes of particles become non-negligible. It is however difficult 
to estimate the lower limit of L, because the finite-volume effect is rather uncontrollable. We shall take the spatial 
extents L = 8, 10,12,16 at /3 = 5.7 as a trial. 

We take periodic boundary conditions in all directions for the gauge field, whereas we impose periodic boundary 
conditions on the spatial directions and the Dirichlet boundary condition on the temporal direction for the quark field 
in order to avoid possible contaminations from those propagating beyond the boundary at t = 0 in (anti)periodic 
boundary conditions. Since the source of possible contaminations is peculiar to the pentaquark and has not been 
properly noticed in previous studies, it is worthwhile to dwell on this problem for a moment. 
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Let us denote the correlators in the pentaquark channel with periodic/antiperiodic boundary conditions and Dirich- 
let boundary conditions as Cp/Ap{T) and Cd{T), respectively. Inserting complete set of states these correlators read 

Cp/Ap{T) = 

m,n 

Cd{T) = (^l«/)(n/|0|m)(m|0|n,)K|D)e-^"^-^"/('^‘-(^+‘"“»-^'‘^‘"=, (10) 

m,nf .rii 

where the states |m), \n), \ni), \nf) are the eigenstates with energies Em, En, E^, En^ respectively. \D) is the 
state which corresponds to the Dirichlet boundary condition and (—)'^” is the factor which represents the ± sign 
with antiperiodic boundary condition. The factor is equal to 1(—1) when \n) contains an even(odd) number 

of valence quarks. It should be noted that the up, down, and strange quark numbers U,D,S for the state \D) are 
restricted to zero since all the quark fields 0) ^'^e set to zero in Dirichlet boundary condition, while there is no 
restriction in the quark sector for the states \n) which appear in periodic/antiperiodic boundary conditions. This 
means that the following states can contribute in the correlators 

\n) = any states including |0), \K), |iV), \NK) ■ ■ ■ for periodic/antiperiodic b.c., 

\nij) = any states(with U = D = S = Q) including |0) • • • for Dirichlet b.c., (11) 

where \K) and \N) are the antiparticle states of \K) and |A^), respectively. Therefore, the correlators have the 
contributions in the long range limit, 

Cp/Ap{T) ~ {Q\Q\NK){NK\Q\{i)e-^^^'^ -f (O|0|0)(0|0|O)e-^®'^ -f (K|0|iV)(iV|0|K)e-®^'^-'®^f("*-^) 

±(iV|0|iL)(ii:|0|iV)e--®«'^--®^^(^*"'^^ ± (lVK|0|O)(O|0|iViL)e--®'^«(^*-^) 

-I- and excited states , (12) 

Cd{T) - {D\Q){Q\Q\NK){NK\Q\Q){0\D)e-^^^'^ + (D|O)(O|0|0)(0|0|O)(O|D)e-®^‘5^ 

-I- and excited states . (13) 

In Fig.^ we give a schematic picture of the contributions to the correlators. The hrst two terms in Eq. m are the 
contributions from the five quark states as given in diagram (A). The third, the fourth and the fifth terms in Eq. 1121 ) 
which correspond to diagrams (C),(D) and (B) are hadronic contributions which propagate beyond the boundary. As 
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FIG. 1: Schematic figure for the explanation on the possible contaminations of the particles propagating over the temporal 
boundary. H{t) and 0(t) are the interpolating operators and their arguments are the distances from the source point. H{t) 
is a generic hadronic operator, which creates and annihilates the particles which cannot decay in quenched QCD. The wave 
lines represent the propagations of states. Resonance states like 0”'" are represented by “5Q” in the figure, as well as NK 
scattering states. The five-quark state can dissociate into forward-propagating nucleon (Kaon) and backward-propagating Kaon 
(nucleon). 

a result, the correlation (0(T -b tsrc)0(tsrc)) inevitably contains unwanted contributions such as 

{K\Q{T + E,,)\N){N\Q{U,,)\K) - (14) 

In this case, the effective mass plot approaches Epi — Ej^ below the NK threshold as T is increased. On the other 
hand, the contributions corresponding to diagrams (C),(D) and (B) do not exist with Dirichlet boundary conditions 
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(Eq.®). Therefore, we find that it would be safest to impose the Dirichlet boundary condition on the temporal 
direction, since no quark can go over the boundary on t=0 in the temporal direction. Although the boundary is 
transparent for the particles composed only by gluons; i.e. glueballs, due to the periodicity of the gauge action, it 
would be however safe to neglect these gluonic particles going beyond the boundary since these particles are rather 
heavy. Then, the correlation (0(T + tsrc)0(isrc)) mainly contains only such terms ((A) in Fig. ^ as 

(vac|0(T + 4rc)|5Q)(5Q|0(<src)|vac) = (1 - 70 ) ^ + (1 + 70 ) ^ ^ (15) 

i i 


with the weight factor and the eigenenergy of i-th state in positive/negative parity channel, respectively. 
One sees that one can now apply the prescription mentioned in the last section. 

One may wonder if these contaminations can be discarded with the parity projection of the correlators by taking 
linear combinations with periodic and antiperiodic boundary conditions. This method indeed works for ordinary three 
quark states where one can single out one of the two contributions diagram (X) and (Y) in Fig.^ However even if 
one takes such linear combinations, one cannot make the contributions from diagram (C) seen in Fa. 11211 cancel out 
as opposed to the contributions from diagram (B) and (D). It is because of the fact that the factor (—for the 
contribution (C) is always equal to 1. (We note here that we can avoid these contaminations using the “averaged 
quark propagator” 113 ■) Some of the previous lattice QCD studies on 0+ adopted a parity projection method using 
the combination with periodic and antiperiodic boundary conditions [T^ 1^ . We stress that one should in principle 
be careful whether the result is free from the contamination owing to the boundary condition which is peculiar to the 
pentaquark and can mimic a fake plateau in the propagator. 

After obtaining the energy spectrum, we carry out a study of the spectral weight for Ks)=( 0.1600,0.1600). 
Introducing two smeared operators 0^smear, 0^smear we compute the following correlators 

C"(T) = ^(0^(X, T + <src)^smear(isrc)), (16) 

X 

from which we extract the spectral weights using a constrained double exponential fit. The details will be explained 
in Sec. IVIII 


IV. LATTICE QCD DATA 

Before obtaining the energies of the lowest state and the 2nd-lowest state, there are only a few simple steps. First, 
we calculate the 2x2 correlation matrix C(T) defined in Eq.@ and obtain the “energies” {Ei(T)} as the logarithm 
of eigenvalues of the matrix product C~^(T + l)C(T). After finding the T range (Tmin <T< Tmax), where 

{Ei{T)} are stable against T, we can extract the energies Ei by the least y-squared fit of the data as Ei = Ei{T) in 
T ■ < T < T 

mm _ _ -^max' 

Since the volume dependence of the energy is crucial to judge whether the state is a resonance or not, a great care 
must be paid in extracting the energy. Therefore, the systematic error by the contaminations from higher excited- 
states should be avoided by a careful choice of fitting ranges. For this reason, we impose the following criteria for the 
reliable extraction of the energies. Although this set of the criteria is nothing more than just one possible choice, we 
believe it is important to impose some concrete criteria for the fit so that we can reduce the human bias for the fit, 
though not completely. 

1. The effective mass plot should have “plateau” for both the lowest and the 2nd-lowest states simultaneously in a 

fit range [Tmin, Tmax], where the length Nfn = — Tmin -I- 1 should be larger than or equal to 3 ( Agt > 3 ). 

2. In the plateau region, the signal for the lowest and the 2nd-lowest states should be distinguishable, so that the 
gap between the central values of the lowest and the 2 nd-lowest energies should be larger than their errors. 

3. The fitted energies should be stable against the choice of the fit range; i.e. the results of the fit with Vgt time 
slices and with Vfit — I time slices should be consistent within statistical errors for both the lowest and the 
2 nd-lowest states. 

4. The lowest state energy obtained by the diagonalization method using the 2x2 correlation matrix should be 
consistent with the value from a single exponential fit for a sufficiently large t. 

If the fit does not satisfy the above conditions, we discard the result since either the data in the fit range may be 
contaminated by higher excited-states or the 2 nd-lowest-state signal is too noisy for a reliable fit. 
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FIG. 2: The “effective mass” plot Ei{T) as the function of T, the separation between source operators and sink operators, in 
= (0, l/2“) channel with the hopping parameters (Ku,d, Ks) = (0.1600.0.1600) on8®x24, 10^x24, 12^x24,16^x24 lattice 
at /3 = 5.7. The stability of each Ei{T) against T means the smallness of the unwanted higher excited-state contaminations. 
The solid line and the dashed lines represent the central value and the error of the fitted masses E^ and E^. 


FigsElshow the “effective mass” plot Ei{T) for the heaviest combination of quarks {Ku,d, Ks) = (0.1600.0.1600). As 
is mentioned in Sec. QH we need to find the T region <T< Tmax) where each Ei{T) shows a plateau. In the 

case of 12^ x 24 lattice in Fig. EJ for example, we choose the fit range of Tmin = 6 and Tmax = 9 {Nf^t = 4). Notice 
that the source operators are put on the time slice with t = tgrc = 6 . The plateau in this region satisfies the above 
criteria so that we consider the fit E^ and E^ for the range 6 < T < 9 as being reliable. The situation is similar for 
the cases of 10^ x 24 and 16^ x 24 lattices. On the other hand, in the case of 8 ^ x 24 lattice we do not find a plateau 
region satisfying the above criteria. 

Figs. shows the “effective mass” plots for the combinations with smaller quark masses. We find that the 

signal is noisier for the lighter quarks and the fit with the smaller volumes 8 ^ x 24 and 10^ x 24 lattices do not satisfy 
the criteria. 


V. LOWEST-STATE ENERGY IN (/, J^) = (0, § ) CHANNEL 
A. The volume dependence of the lowest-state energy 

Now we show the lattice QCD results of the lowest state in / = 0 and ^ channel. The filled circles in Fig.[7| 

show the lowest-state energies Eq in / = 0 and ^ channel on four different volumes. Here the horizontal axis 
denotes the lattice extent L in the lattice unit and the vertical axis is the energy of the state. The lower line denotes 
the simple sum M^ + Mk of the nucleon mass Mn and Kaon mass Mk obtained with the largest lattice. Though 
and Mk are slightly affected by finite volume effects, the deviation of Mn + Mk {L = 8) from Mk + Mk {L = 16) 
is about a few % (Tabled)- We therefore simply use Mk + Mk {L = 16) as a guideline. 

At a glance, we find that the energy of this state takes almost constant value against the volume variation and 
coincides with the simple sum Mk + Mk- We can therefore conclude that the lowest state in / = 0 and ^ 

channel is the NK scattering state with the relative momentum |p| = 0. The good agreement with the sum Mk + Mk 
implies the weakness of the interaction between N and K. In fact, the scattering length in the 1 = 0 channel is known 
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FIG. 3: The “effective mass” plot Ei{T) as the function of T, the separation between source operators and sink operators, in 
= (0, l/2“) channel with the hopping parameters {ku, d,Ks) = (0.1600.0.1650) on 8® x 24, 10® x 24, 12® x 24,16® x 24 
lattice at f3 = 5.7. 
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FIG. 4: The “effective mass” plot Ei{T) as the function of T, the separation between source operators and sink operators, in 
= (0,1/2”) channel with the hopping parameters {ku, d,K,s) = (0.1650.0.1600) on 8® x 24, 10® x 24, 12® x 24,16® x 24 
lattice at f3 = 5.7. 
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FIG. 5: The “effective mass” plot Ei{T) as the function of T, the separation between source operators and sink operators, in 
= (0, l/2“) channel with the hopping parameters {ku, d,Ks) = (0.1625.0.1650) on 8® x 24, 10® x 24, 12® x 24,16® x 24 
lattice at f3 = 5.7. 
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FIG. 6: The “effective mass” plot Ei{T) as the function of T, the separation between source operators and sink operators, in 
= (0,1/2”) channel with the hopping parameters {ku, d,K,s) = (0.1650.0.1650) on 8® x 24, 10® x 24, 12® x 24,16® x 24 
lattice at f3 = 5.7. 
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to be tiny ( {I = 0) = —0.0075 fm ) from compilations of hadron scattering experiments [s^ whereas the 
current algebra prediction from PC AC with SU(3) symmetry predicts that the scattering length {I = 0) = 0. 


B. Comparison with the previous lattice work 

We here compare our data with the previous lattice QCD studies, which were performed with almost the same 
conditions as ours, in order to confirm the reliability of our data. 

The well-known hadron masses mT^^rriK and mjv listed in TableOlcan be compared with the values in Ref. Our 
data are consistent with those in Ref. The lowest NK scattering state in (/, J^) = (0, ^ ) channel is carefully 
investigated in Ref. with almost the same parameters as our present study. It is worth comparing our data 
with them. For the complete check of our data, we re-extract the lowest state by the ordinary single-exponential fit 
of the correlator as (0(T -|- tsrc)0waii(fsrc)) = in the large-T region, and compare them with the present 

lattice data Eq obtained by the multi-exponential method as well as the data in Ref. . In Table Ql we list 
the data of the lowest state E^k obtained by the single-exponential fit. They almost coincide with the present 
data Eq extracted by the multi-exponential method with about 1% deviations, which may be considered as the 
slightly remaining contaminations of the higher excited-state. In Ref. , the authors extracted the energy difference 
SE = Enk — {En + Ek) with the hopping parameters Ku,d,s = 0.1640 using 12^ x 20 lattice. We therefore compare 
our data SE = Enk — {En + Ek) obtained with the hopping parameters {nu,d^i^s) = (0.1650,0.1650) on 12^ x 24 
lattice. The energy difference SE in our study is found to be SE = —0.0128(38), which is consistent with the value of 
—0.0051(38) in Ref. pSi taking into account that this error includes only statistical one. 

It is now confirmed that the lowest state extracted using the multi-exponential method is consistent with the 
previous works and that our data and method are reliable enough to investigate the 2nd-lowest state in this channel. 


VI. 2ND-LOWEST STATE ENERGY IN (1, J^) = (0, i“) CHANNEL 

The (J, J^) = (0, ^ ) state is one of the candidates for 0+(154O). Since 0+ is located above the NK threshold, it 
would appear as an excited state in this channel. We show the lattice data of the 2nd-lowest state in this channel. 

In order to distinguish a possible resonance state from NK scattering states, we investigate the volume dependence 
of both the energy and the spectral weight of each state. It is expected that the energies of resonance states have 

small volume dependence, while the energies of NK scattering states are expected to scale as ^Mfj + |^np -|- 

\J+ l^fip according to the relative momentum between N and K on a finite periodic lattice, provided that 
the NK interaction is weak and negligible. We can take advantage of the above difference for the discrimination. 


A. Possible corrections to the volume dependence of NK scattering states 

A possible candidate for the volume dependence of the energies of NK scattering states is the simple formula as 
E^k{L) = -|- l^np 4- M|. -|- Iwith the relative momentum between N and K in finite periodic 

lattices, which is justified on the assumption that nucleon and Kaon are point particles and that the interaction 
between them is negligible. In practice, there may be some corrections to the volume dependence of E^j^{L). We 
therefore estimate here three possible corrections; the existence of the NK interaction, the application of the momenta 
on a finite discretized lattice and the estimation of the implicit finite-size effects. 

There can be small hadronic interactions between nucleon and Kaon, which may lead to correction to naively 
expected energy spectrum E^i^{L) of the NK scattering states. Using Liischer formula [s^, one can relate the 
scattering phase shift to the energy shift from iJ”^(L) on finite lattices. For example, in the case when a system 
belongs to the representation Af of cubic groups, which is relevant in the present case, the relation between the phase 
shift and the possible momentum spectra is 


g2z5o(fe) ^ 2{l;q^)+iTTlq 
Z(l; q^) — iiriq 


( 17 ) 
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Here Z{s] q^) is the Zeta function defined as 


Z{s\q^) 


1 


•\/4^ 


E 


1 

(n^ — ’ 


(18) 


with the eigenenergy g on a finite lattice. We have simply omitted the corrections from the partial waves with angular 
momenta higher than the next smallest one (1 = 4). Although our current quark masses are heavier than those of 
the real quarks, we use the empirical values of the phase shift in NK scattering in Ref.j^, by simply neglecting the 
quark mass dependence. The correction using the empirical values results in at most a few % larger energy than the 
simple formula within the volume range under consideration; the energies are slightly increased by the weak 

repulsive force between nucleon and Kaon. 

One may claim that one has to adopt momenta on a finite discretized lattice: |pp = 4sin^(^) • |np for Kaon and 
IpP = sin^(^) • |np for nucleon, respectively. This correction turns out to be within only a few % lower energy than 
although it is not certain whether this correction is meaningful or not for composite particles like nucleon 

or Kaon. 

We find that these corrections lead to at most a few % deviations from We then neglect these corrections 

for simplicity in the following discussion and use the simple form E^j^{L). 


So far, we have neglected the implicit finite-size effects in Efjj^{L), other than the explicit ones due to the lattice 

momenta |p| = ^|n|. Some smart readers may suspect that the dispersions Mfj + |^np and \J 

may be affected by the uncontrollable finite-size effects due to the finite sizes of N and K, and no longer valid. In 
order to make sure of the small implicit artifacts especially with |n| = 1, which we are mainly interested in, we also 


calculate the sum ^ -|- ^ of energies of nucleon a]"* ^ 






and Kaon E 


|n| = l 
K 


momentum |p| = ^|n| = 2t: jL. We extract E^ and E^ from the correlators 


with the smallest non-zero lattice 
‘T«'*(iV(x,t-|-tsrc)A^(0,tsrc)) and 


‘^{K{5c, t+tsic)K{0, tsrc)). These results are denoted by the open squares in Fig. das the sum E',^' ‘'EE 


The upper lines in Fig. dshow ^ -I- y/Mj. + {2 'k/LY. The deviations of ^ + E 


src ) 

|n| = l ptn| = l 

N +F'_fs- 
,|n| = l 


p|nl = l 

from are very small, which implies the smallness of the implicit finite-size artifacts. Therefore, provided that 

the interaction between nucleon and Kaon is weak, which we assume throughout the present analysis, the naive 
expectation for the 2nd-lowest NK scattering states denoted by the upper line in Fig. d would be able to follow the 
energies of the 2nd-lowest NK scattering states even on the L = 8 lattices in our setup. 


B. The volume dependence of the 2nd-lowest state energy 

We compare the lattice data E^ with the expected behaviors for the 2nd-lowest NK scattering states. 

The filled squares in Fig. d denote E^, the 2nd-lowest-state energies in this channel. The black and gray symbols 
are the lattice data obtained by the fits with Nat — 1 and A^fit time slices, respectively (see the criterion.3 in the 
Sec.lTV|l. The upper line shows the expected energy-dependence on V of the 2nd-lowest NK scattering state ^ = 
\/+ (27r/F)2 -I- {2ttILY estimated with the next-smallest relative momentum between N and K, and 

with the masses Mk and Mm extracted on the F = 16 lattices. Although the lattice QCD data E^ and the expected 
lines almost coincide with each other on the F = 16 lattices, which one may take as the characteristics of the 

2nd-lowest scattering state, the data E^ do not follow in the smaller lattices. (At the smallest lattices with 

F = 8 (1.4 fm) in the physical unit, some results apparently coincide with each other again. However we consider 
that the volume with F ~ 1.4 fm is too small for the pentaquarks; it is difficult to tell which is the origin of the 
coincidence, uncontrollable finite volume effects of the pentaquarks or expected volume dependence of the 2nd-lowest 
NK scattering state.) Especially when the quarks are heavy, composite particles will be rather compact and we expect 
smaller finite volume effects besides those arising from the lattice momenta p = ^n. Moreover the statistical errors 
are also well controlled for the heavy quarks. Thus, the significant deviations in 1.5 F ^ 3 fm with the combination 
of the heavy quarks, such as (nu,d, Ks) = (0.1600,0.1600), are reliable and the obtained states are difficult to explain 
as the NK scattering states. 

Therefore one can understand this behavior with the view that this state is a resonance state rather than a scattering 
state. In fact, while the data with the lighter quarks have rather strong volume dependences which can be considered 
to arise due to the finite size of a resonance state, the lattice data exhibit almost no volume dependence with the 
combination of the heavy quarks especially in 1.5 ^ F < 3 fm, which can be regarded as the characteristic of resonance 
states. 
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FIG. 7: The black (gray) filled-squares denote the lattice QCD data of the 2nd-lowest state in (I, J^) = (0, | ) channel 
extracted with — 1 data plotted against the lattice extent L. The Hlled circles represent the lattice QCD data of 

the lowest state in (7, J^) = (0, ^ ) channel. The open symbols are the sum of energies of nucleon and 

Kaon 75^'^^ with the smallest lattice momentum |p| = ^|n| = 27r/L. The upper line represents r]p[2 + Mf, + |pP 

with IpI = 27r/L the smallest relative momentum on the lattice. The lower line represent the simple sum Mm + Mk of the 
masses of nucleon Mm and Kaon Mk- We adopt the central values oi Mm and Mk obtained on the largest lattice to draw the 
two lines. 
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VII. THE VOLUME DEPENDENCE OF THE SPECTRAL WEIGHT 

For further confirmation, we investigate the volume dependence of the spectral weight [^. As mentioned in Sec. nn 
the correlation function (O(T)Ol(O)) can be expanded as (O(r)Ol(O)) = ^ WiC . The spectral weight of the i-th 
state is defined as the coefficient Wi corresponding to the overlap of the operator 0(t) with the i-th excited-state. 
The normalization conditions of the field ip and the states \i) give rise to the volume dependence of the weight factors 
Wi in accordance with the types of the operators 0{t). 

For example, in the case when a correlation function is constructed from a point-source and a zero-momentum 
point-sink, as + isrc)0(O,isrc)), the weight factor Wi takes an almost constant value if |i) is the resonance 

state where the wave function is localized. If the state |i) is a two-particle state, the situation is more complicated. 
Nevertheless if there is almost no interaction between the two particles, the weight factor is expected to be proportional 
to y. In the case when a source is a wall operator ©waiif^src) as taken in this work, a definite volume dependence 

of Wi is not known. Therefore, we re-examine the lowest state and the 2nd-lowest state in (J, J^) = (0, i ) channel 
using the locally-smeared source 0^smear(isrc) = Sser tsrc) with F = ({0, 3}, {0, 3}, {0, 3}), which we introduce 
to partially enhance the ground-state overlap. Since smeared operators, whose typical sizes are much smaller than the 
total volume, can be regarded as local operators, we can discriminate the states using the locally-smeared operators 
as in the case of point operators. (We also investigated the weight factor using the point source. The results are 
consistent with those obtained using the locally-smeared one, but are rather noisy.) We adopt the hopping parameter 
{Ku,d,Ks) = (0.1600,0.1600) and additionally employ 14^ x 24 lattice for this aim. 




Volume V [lattice unit] Volume V [lattice unit] 

FIG. 8: The spectral weight factors defined in Sec. IVTTI are plotted against the lattice volume V. The left figure shows Wo 
for the lowest state in (/, J^) = (0, | ) channel and the right figure shows Wi for the 2nd-lowest state. In the case when the 
weight factor Wi for the i-th state |i) in a point-point correlator shows no volume dependence, \i) is likely to be a resonance 
state. On the contrary, when the i-th state |i) is a two-particle state, Wi scales according to 1/V. 

We extract Wq and Wi using the two-exponential fit as + 4rc)©^smear(0, ^src)) = Woe~^°'^ + Wie~'^^'^. 

The fit with four free parameters Wq, Wi, Vq and Vi is however unstable and therefore we fix the exponents using 
the obtained values Eq and . The weight factors Wq and Wi are then obtained through the two-parameter fit as 
^-(©^(x, T + tsrc)©^smear(0, tgrc)) = Woe~^o ^ -|- VFie“^i ^ in as large t range (T^in < T -|- Urc < Tmax) as possible 
in order to avoid the contaminations of the higher excited-states than the 2nd-excited state (3rd-lowest state), which 
will bring about the instability of the fitted parameters. The fluctuations of Eq and E^ are taken into account 
through the Jackknife error estimation. Fig. |H| includes all the results with the various fit range as (Tniin,Tniax)= 
(16,18),(17,19),(18,20) to see the fit-range dependence. Though the results have some fit-range dependences, the 
global behaviors are almost the same among the three. 

The left figure in Fig. [S] shows the weight factor Wq of the lowest state in (/, J^) = (0, ^ ) channel against the 
lattice volume V. We find that Wq decreases as V increases and that the dependence on V is consistent with A, 
which is expected in the case of two-particle states. It is again confirmed that the lowest state in this channel is the 
NK scattering state with the relative momentum |p| = 0. Next, we plot the weight factor Wi of the 2nd lowest state 
in the right figure. In this figure, almost no volume dependence against V is found, which is the characteristic of 
the state in which the relative wave function is localized. (In Appendix, we try another prescription to estimate the 
volume dependences of the spectral weights, which requires no multi-exponential fit.) This result can be considered 
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as one of the evidences of a resonance state lying slightly above the NK threshold. 

To summarize this section, the volume dependence analysis of the eigenenergies and the weight factors of the 2nd- 
lowest state in = (0, i ) channel suggests the existence of a resonance state. Although there remain the 

statistical errors and the possible finite-volume artifacts, the data can be consistently accounted for assuming the 
2nd-lowest state to be different from ordinary scattering states. If the 2nd-lowest state were an ordinary scattering 
state, one had to assume a large systematic errors for heavier quarks which is hard to understand consistently. 


VIII. (J, J^) = (0, CHANNEL 






Extent [lattice unit] 


Extent [lattice unit] 


FIG. 9: The lattice QCD data in the (I, J^) = (0, channel are plotted against the lattice extent L. The solid line denotes 
the simple sum Mjv* -I- Mk of the masses of the lowest-state negative-parity nucleon Mn* and Kaon Mk obtained with the 
largest lattice. 


In the same way as (/, J^) = (0, ^ ) channel, we have attempted to diagonalize the correlation matrix in (/, J^) = 

( 0 , channel using the wall-sources 0waii(i) and the zero-momentum point-sinks ^^is channel, the 

diagonalization is rather unstable and we find only one state except for tiny contributions of possible other states. 
We plot the lattice data in Fig. 0 One finds that they have almost no volume dependence and that they coincide 
with the solid line which represents the simple sum Mm* + Mk of Mm* and Mk, with Mm* the mass of the ground 
state of the negative-parity nucleon. From this fact, the state we observe is concluded to be the N*-K scattering state 
with the relative momentum |p| = 0. It may sound strange because the p-wave state of N and K with the relative 
momentum |p| = 27r/L should be lighter than the N*-K scattering state with the relative momentum |p| = 0; this 
lighter state is missing in our analysis. This failure would be due to the wall-like operator 0waii(i)- The fact that the 
wall operator 0 waii(i) is constructed by the spatially spread quark fields with zero momentum may leads to 

the large overlaps with the NK scattering state with zero relative momentum. The relation between operators and 
overlap coefficients is an interesting problem and is to be explored in detail for further studies. Anyway, the strong 
dependence on the choice of operators suggests that it is needed to try various types of operators before giving the 
final conclusion. 

Before closing this section, we show the spectral weight ITq^ obtained by the fit using the form (0^(t -I- 
4rc)0^waii(4rc)) = W^e~^o* in Fig. ^1 Although one sees the -^-like volume dependence in Fig. ^1 one can 
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Volume V [lattice unit] 


FIG. 10: The spectral weight factor of the extracted state in (1, J^) = (0, channel with the hopping parameters {Ku,d, Ks) = 
(0.1600, 0.1600) is plotted against the lattice volume V. We note here that the l/l/-like volume dependence is not always the 
characteristics of scattering states when we do not adopt point-point correlators as is the present case. To conclude from this 
dependence, we need to determine the precise volume dependences of the weight factors in wall-point correlators. 


conclude nothing only from this behavior unless the precise volume dependences of the weight factors in wall-point 
correlators are estimated. 


IX. DISCUSSION 
A. Operator dependence 

We here mention the operator dependences in (I, J^) = (0, i ) channel. As is seen in Sec. IVIIII the overlap 
factors with states strongly depend on the choice of operators. We survey the effective masses of the five correlators; 

( Q ^ (x, r + 4rc ) Q Vail (4rc ) ) , V V + 4rc)0^wall (Vc)) , Z]x(0Hx, ^-l-tsrc)© V Vc)) , Z]x(0^(x, T 

and X]x(0^(^’ ^ + 4rc)0^(4rc))- The effective mass E(T) is defined as a ratio between correlators with the temporal 
separation T and T -|- 1, 


E{T) = In 


(O(r)O(0)t) 
(o(r + i)o(o)t)’ 


which can be expressed in terms of the eigenenergies and spectral weights as 


E{T) = In 


E, lU,e-sdT+i) 


Eq + 


(E^-Eo)T , 

W 2 


(19) 


( 20 ) 


A plateau in E(T) at Eq implies the ground-state dominance in the correlator. Effective mass plots E{T) are often 
used to find the range where correlators show a single-exponential behavior; the higher excited-state contributions 
{i > 0) are negligible in comparison with the ground-state component 
Here, 0^ is an interpolation operator defined as 

03(x) = V‘’V“fe*’8*'[ne(x)Cd/(x)]K(x)C75V(^)]t^V(x) (21) 

which has a di-quark structure similar to that proposed by Jaffe and Wilczek |^, and is also used in Refs. [T^l^l^ . 
Eig.^^Qshows the effective mass plots constructed from ^ )—(0 (x, T-t-tsi.(;)0 waii(^src))) ^ (0 (x, T-t-tgj.(,)0 wall (^src)) ^ 

+ tsrc)O^(tsrc)}, + tsrc)0^(Vc)) and J^ji(0^(x,T + tsrc)O^(tsrc)}- One can see two typical 

behaviors in the figure. One is the line damping from a large value to the energy E^^ of the lowest NK scattering 
state. The other is the one arising upward to Eq . Surprisingly, the differences of the spinor structure or the color 
structure among the operators are hardly reflected in the effective mass plots. The difference is enough to perform 
the variational method but seems insufficient for a clear change of the effective mass plots. Instead, the effective mass 
plots seem sensitive to the spatial distribution of operators. The upper three symbols are data using the point-point 
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FIG. 11: The effective mass plots constructed from ^j(0^(x,T + tsrc)0^waii(tsrc)), ^j(0^(x,r + tsrc)0^waii(tsrc)), 

^j(0^(x,T + tsrc)0H^src)), + tsrc)0^(tsrc)) and ^j(0®(x,T + tsrc)0^(tsrc)), in {1,J^) = (0, f”) channel in 

16^ X 24 lattice at /3 = 5.7 employing the hopping parameters «:s)=(0.1600,0.1600) are plotted, along with the dashed 

line which denotes the lowest-state energy . 


correlators and the lower two symbols are those from the wall-point correlators. This means the overlap factor with 
each state is controlled mainly by the spatial distribution rather than the internal structure of operators, except for 
the overall constant. The spatially smeared operators seem to have larger overlaps with the scattering state with the 
relative momentum |p| = 0. (One can find that the overlap factor of the wall operator with the observed state in 
(/, J^) = (0,1 ) in Fig.^lis 1000 times larger than those of point operators in Fig.|Sl) One often expects that the 
overlap with a state could be enhanced using an operator whose spinor or color structures resembles the state. We 
find however no such tendency in the present analysis. The insensitivity to the spinor structures may come from the 
fact that the KN-type operator (0^) and the di-quark type operator (0^) are directly related by a factor of 75 and a 
Fierz rearrangement j22l |. Though we have no idea about the mechanism of the insensitivity to the color structure at 
present, this insensitivity would have some connection with the internal color-structure of 0 '*". 

The upper three data slowly damp and do not reach the lowest energy Eq in this T range, which can be explained 
in terms of the spectral weight. As is seen in Fig. O Wq is ten-times smaller than Wi in the case of the point-point 
correlator. Then, the term in the effective mass survives at relatively large T. Hence the effective 

mass needs larger T to show a plateau at Eq . The insensitivity of the overlaps to the internal structure of operators 
could be helpful for us: We have adopted two operators whose color and spinor structures are different from each 
other. Although the difference is enough in (/, J'^) = (0, ^ ) channel, it may be insufficient in (/, J^) = (0, 
channel, which leads to the failure in the diagonalization. If we use operators with spatial distributions different from 
each other, it would be more effective in the diagonalization method. 


B. Comment on other works 

Here we comment on other works previously published, especially for the pioneering works by Csikor et al. |1S| and 
Sasaki |19|. The simulation condition for the former is rather similar to ours. 

Csikor et al. first re por ted the possible pentaquark state slightly above the NK threshold in (/, J^) = (0, ^ ) 
channel in 0 In Ref. |l 8 | , they tried chiral extrapolations and taking the continuum limit at the quenched level for 
the possible pentaquark state. However they used the single-exponential fit analysis for the non-lowest state, namely 
the possible pentaquark state, for the main results. It is difficult to justify their result unless the coupling of the 
operators to the lowest NK state is extremely small. 

Sasaki found a double-plateau in the effective mass plot and identified the 2nd-lowest plateau as the signal of 0+. 
Unfortunately, we does not find a double-plateau in the present analysis. The double-plateau-like behavior in effective 
mass plots can appear only under the extreme condition that Wi is much larger than Wq. Wi which is ten times 
larger than Wq seen in Fig. 0 and the statistical fluctuations may cause the deviation of the effective mass plot from 
the single monotonous line. In fact, the effective mass plot very slowly approaches Eq as T increases in Fig. He 
extracted the mass of the next-lowest state with a single and double exponential fits. The result do not contradict 
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with ours. 

Ref. 1^ reports a lattice QCD study which adopted the overlap-fermions with the exact chiral symmetry. The 
hybrid-boundary method was suggested in Ref. and the authors tried to single out the possible resonance state. 
In these two studies, the absence of resonance states with a mass a few hundred MeV above the NK threshold was 
concluded. We have not found the resonance state which coincides just with the mass of 0'*' in the chiral limit. In 
this sense, the results in Refs. |23,|23 are not inconsistent with ours. 


C. chiral extrapolation 

We perform chiral extrapolations for Kaon, nucleon, NK threshold (a simple sum of a Kaon mass and a nucleon 
mass) and the 2nd-lowest state in the (/, J^) = (0, i ) channel. We adopt the lattice data with 16^ x 24 lattice, the 
largest lattice in our analysis. One can find in Fig. |3 that the 2nd-lowest state, which is expected to be a resonance 
state, is already affected by the finite volume effects for L < 12 with the lightest combination of quarks, and we 
therefore adopt the largest-lattice data for safety. We can expect from this fact that the typical diameter of this 
resonance is about 2 fm or longer and that it is desirable to use larger lattices than, (2.5 fm)^ for the analysis of 0+. 

In Fig. [El Mk + Mm and obtained with each combination of quark masses for 12^ x 24 and 16^ x 24 lattices 
are plotted against mj. We assume the linear function of quark masses, = 6oo + &io''7ru,d + ^oi^sj for 

nucleon and the 2nd-lowest state with 6^ free parameters fitted using the five lattice data. We determine the critical k 
(kc) by and fix the Kg so that the physical Kaon mass is reproduced in the chiral limit, using the form for pseudo 
scalar mesons Epg(rnu,d, nTs) = aiQUiu^d + ooiWg. The chiral-extrapolated values of Mk, Mm, Mk + Mm and E^ 
for 16^ X 24 lattice are 0.4274(12), 0.7996(60), 1.227(6) and 1.500(52) in the lattice unit and 0.5001(14), 0.9355(70), 
1.436(7) and 1.755(61) in the unit of GeV, respectively. We find that the results for 12^ x 24 lattice are consistent 
within errors as shown in Fig. 1121 

The value of Kj“=1755(61) MeV in the chiral limit is significantly larger than the mass of 0+(154O) in the real 
world. How can we interpret this deviation? One possibility is the systematic errors from the discretization, the 
chiral extrapolation, or quenching. Another possibility is that the observed 2nd-lowest state might be a signal of a 
resonance state lying higher than O'*'. Unfortunately there is no clear explanation at this point. Obviously more 
extensive studies on finer lattices with lighter quark masses in unquenched QCD are required. However, we can at 
least conclude that our quenched lattice calculations suggest the existence of a resonance-like state slightly above the 
NK threshold for the parameter region we have investigated. 



(mj^ [lattice unit] (mj^ [lattice unit] 


FIG. 12: A comparison of the chiral extrapolations in (7, J^) = (0, | ) channel on 12® x 24 (left) and 16® x 24 (right) lattices 
for the NK threshold energy Mk + Mm and the 2nd-lowest state energy . Mk + Mn and are plotted against M^. 
The filled circles (triangles) denote the energies of the NK threshold (2nd-lowest state) in the chiral limit. Ks is fixed to be 
Ka ~ 0.1620 so that the physical Kaon mass is reproduced in the chiral limit. 
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X. SUMMARY AND FUTURE WORKS 

We have performed the lattice QCD study of the {S,I,J) = (+1,0, states on 8^ x 24, 10^ x 24, 12^ x 24 and 
16^ X 24 lattices at /3=5.7 at the quenched level with the standard plaquette gauge action and Wilson quark action. 
To avoid the possible contaminations originating from the (anti)periodic boundary condition, which are peculiar to 
the pentaquark and have not been properly noticed in previous studies, we have adopted the Dirichlet boundary 
condition in the temporal direction for the quark field. With the aim to separate states clearly, we have adopted two 
independent operators with 1 = 0 and ^ so that we can construct a 2 x 2 correlation matrix. 

From the correlation matrix of the operators, we have successfully obtained the energies of the lowest state and the 
2nd-lowest state in the (/, J^) = (0, i ) channel. The volume dependence of the energies and spectral weight factors 
show that the 2nd-lowest state in this channel is likely to be a resonance state located slightly above the NK threshold 
and that the lowest state is the NK scattering state with the relative momentum |p| = 0. As for the (/, J^) = (0, 
channel, we have observed only one state in the present analysis, which is likely to be a N*K scattering state of the 
ground state of the negative-parity nucleon N* and Kaon with the relative momentum |p| = 0. 

We have also investigated the overlaps using five independent operators. As a result, we have found that the 
overlaps seem to be insensitive to the spinor and color structure of operators while the overlaps are mainly controlled 
by the spatial distributions of operators, at least for a few low-lying state in this analysis. For the diagonalization 
method, it may be more effective to vary the spatial distributions rather than the internal structures. 

The volume dependence of Af suggests that this resonance-like state in the (/, J^) = (0, ^ ) channel is a rather 

spread object with the radius of about 1 fm or more. The possibility of a resonance state lying in (/, J^) = (0, ^ ) 
channel is desired to be confirmed by other theoretical studies, such as quark models, QCD sum rules, string models 
and so on ESIHIIllIllilllllli. Unfortunately, four quarks uudd and one antiquark s in J = ^ state can hardly 
reproduce the unusually narrow width of 0'*' so far, while the obtained mass in = i channel could be assigned 

to the observed resonance state Hence, = | or states are favored to reproduce the width in 

terms of the quark model. However, there are many unknown problems left so far such as the internal structure of 
multi-quark hadrons 0, |3^ 0) E3 dynamics of the string/flux-tubes The discovery of 0+ 

gives us many challenges in the hadron physics and more detailed theoretical study including the lattice QCD studies 
are awaited. 

For further analyses, a variational analysis using the 3x3 correlation matrix or larger matrices will be desirable. 
The observation of wave functions will be also useful to distinguish a resonance state from scattering states and to 
investigate the internal structures of hadrons. We can use the lattice QCD calculations in order to estimate the decay 
width and to study the flux-tube dynamics which should give useful inputs for model calculations. 
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APPENDIX A: ADDITIONAL ESTIMATIONS OF WEIGHT FACTORS 

In this appendix, we make another trial to estimate volume dependences of weight factors in (J, J^) = (0, ^ ) 
channel. As seen in Sec. ivni we have extracted the weight factors using double-exponential fit, which is however 
rather unstable and we have therefore fixed the exponents. We here discuss the possibility of methods without any 
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_(K^.d,K4=(0.1650,0.1650)_ 

size AAx A/iv Mn* Eq Enk E^ 

8=^ X 24 0.4378(49) 0.4378(49) 1.0463(89) 1.4706(338) — — 1.3841(202) — 

10® X 24 0.4543(17) 0.4543(17) 1.0591(91) 1.4313(292) — — 1.4704(202) 2.0987(45) 

12® X 24 0.4563(13) 0.4563(13) 1.0281(74) 1.4760(309) 1.4601( 75) 1.7881(829) 1.4715(107) 2.0502(26) 

16® X 24 0.4556(11) 0.4556(11) 1.0143(44) 1.4791(278) 1.4616( 56) 1.7157(452) 1.4743( 74) 1.9951(19) 


_(K^,d,K«)=(0.1625,0.1650)_ 

size AItt Afx Afjv A7jv* Eq E-^ E^ 

8® X 24 0.5747(21) 0.5130(31) 1.2112(133) 1.5973(220) 1.6123(139) 2.6447(3322) 1.6119(195) 2.2065(368) 

10® X 24 0.5785(11) 0.5199(14) 1.1814( 55) 1.5712( 93) 1.6673( 96) 2.0912(2305) 1.6687(127) 2.2120(393) 

12® X 24 0.5792(10) 0.5205(10) 1.1655( 47) 1.5930(175) 1.6612( 55) 1.9228(348) 1.6763( 70) 2.2145(293) 

16® X 24 0.5789( 9) 0.5209(10) 1.1590( 37) 1.5888(169) 1.6636( 42) 1.8769(225) 1.6745( 50) 2.1734(223) 


_ (K^.d,Ai:«) = (0.1600,0.1650) _ 

size AItt A/x Afx Mn* Eq E^ E^k E^ 


8® X 24 0.6839(18) 0.5761(25) 1.3270(110) 1.6956( 54) 1.8002(121) 2.5420(1226) 1.8109(259) 2.3666(122) 

10® X 24 0.6873(10) 0.5819(12) 1.3070( 48) 1.6810(131) 1.8549( 78) 2.1797(1067) 1.8574( 90) 2.3547(253) 

12® X 24 0.6883( 9) 0.5816(11) 1.3003( 48) 1.7198(221) 1.8627( 51) 2.0736(306) 1.8680( 57) 2.4004(396) 

16® X 24 0.6867( 9) 0.5818( 9) 1.2921( 30) 1.7193(222) 1.8546( 37) 2.0429(156) 1.8668( 42) 2.3403(162) 


_(K^,dA«)=(0-1650,0.1600)_ 

size Afvr Mk Mn Mn* Eq E^ Enk Eq 

8® X 24 0.4378(49) 0.5761(25) 1.0463( 89) 1.4705(338) — — 1.5150(259) 2.0510(732) 

10® X 24 0.4543(17) 0.5823(17) 1.0774(128) 1.4313(292) — — 1.6088(173) 2.1897(1549) 

12® X 24 0.4563(13) 0.5816(11) 1.0281( 74) 1.4760(309) 1.5838( 72) 1.8368(538) 1.5951(100) 2.1645(440) 

16® X 24 0.4556(11) 0.5818( 9) 1.0143( 44) 1.4791(278) 1.5852( 55) 1.7855(313) 1.5987( 73) 2.0823(176) 


_(K,t,d,Ka) = (0.1600,0.1600)_ 

size Mtt Mk Mn Mn* Eq E^ Enk Eq 

8® X 24 0.6839(18) 0.6839(18) 1.3270(110) 1.6956( 54) — — 1.9239(226) 2.4376(105) 

10® X 24 0.6873(10) 0.6873(10) 1.3070( 48) 1.6810(131) 1.9622( 73) 2.2085(478) 1.9617( 87) 2.4319(215) 

12® X 24 0.6883( 9) 0.6883( 9) 1.2987( 42) 1.7198(221) 1.9632( 51) 2.1528(195) 1.9705( 56) 2.4820(312) 

16® X 24 0.6867( 9) 0.6867( 9) 1.2906( 30) 1.7193(222) 1.9641( 40) 2.1158(153) 1.9723( 41) 2.4180(138) 


TABLE I: The pion masses A7.,r, Kaon masses Mk, nucleon masses Mn, masses of the ground state of negative-parity nucleon 
Mn* , energies of the lowest state Eq , energies of the 2nd-lowest state E^, energies of the lowest state Enk (obtained by single¬ 
exponential fit) in the (7, J^) = (0, ^ ) channel are listed. The energies of the obtained state Eq in the (7, J^) = (0, |^) 
channel are also listed. and Ks are the hopping parameters for u, d quarks and s quark respectively. 


multi-exponential fits. Let us again consider N x N correlation matrices constructed by 7 -sink a-source and 7 -sink 
/3-source correlators. Here a, (3 ,7 denote the types of operators, such as “point” or “wall” or “smear” and so on. The 
notations are the same as those in Sec. im The 7 — a and 7-/3 correlation matrices are described as 

C7“(T) = {Q\{T + 4rc)^(7src)) = (C''>'U(r)C“)/j + + ... (Al) 

(T) = (07(r + tsrc)^(7src)) = + ... (A2) 

with N X N matrices -I- ...) and [d'jje~^^"'^ + ...) being possible higher excited-state contaminations. We 

hereby consider two quantities; \C'lj{T)C'Jj{T + 1)“^]^C77(T) defined using one type of the correlation matrix and 
[C'Jj{T)]~^C'Jj{T), which with large T lead to 

[c 77 (T)c 77 (r -b -b e{d{t)) +.... 


(A3) 
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FIG. 13: The left figure shows det y\Cfj{T)Cfj{T + 1) as the function of T on each volume. The middle 

figure is the plot of det ^ + 1)”^] ^ Cfj'^(r)^ against T. Det ([Cfj^(r)]“^CfJ(r)^ is plotted in the right figure. 

Det([Cf|(T)CF|(r+l)-i]^Cf|(T)), det ([Cf7(r)Cf7(r+l)-i]^Cf7(T)) and det([Cf7(r)]-icF|(r)) show plateaus 
in the large T region and coincide with det(C^^C'®), det(C'^tC'''^) and det((C'^)“^C®) respectively. 


and 


[C]j{T)]-^C]^(T) = + T'{D{T)) + ..., (A4) 

respectively. Here !F{D{T)) and J-'{D{T)) are terms including N x N diagonal matrix D{T) = 
diag(e“(^"“^^“i)^,...,Then, each component of [C]j{T)C]j{T+l)~^]'^C]j{T) and 
[C]j(T)]~^C]j{T) gets stable and shows a plateau in large T region, where T{D{T)) and T'{D{T)) are negli¬ 
gible. 

Next, we relate these quantities to spectral weights. For this aim, we simply take the determinants. In the 
case when the correlation matrices are 2x2 matrices, the determinant det is explicitly written as 

detC^/det (7“ = C^jC^j/e^' ''' Cqj,C^j, , and the determinant det is expressed as (detC^'f) x (detC^) = 

{e^C]lC'j\) X CqjiC^j,) = e^' C'JICqj,C'j\C^j, . The term C^jC^j (CqjC^j) denotes the product of the 

overlaps of the /3(a)-type operator with the lowest state and the 2nd-lowest state. On the other hand, C]qCqj, 
(C'jJC'fj,) corresponds to the spectral weight for the lowest (2nd-lowest) state in the P-j correlator in terms of a 
volume dependence. 

Let us consider the several cases when (a,/3, 7 )={W(wall), S(smeared), P(point)}. The term det(C'^'l'C'®) behaves 
showing the same volume dependence as the product of the spectral weights for the lowest and the 2nd-lowest state in 
the smeared-point correlator, which should be^ y xl = y if the lowest state is a scattering state and the 2nd-lowest 

state is a resonance state. The left panel in Fig. ^Irepresents det ^[CFj (T)CFj (T -|- 1)“^]^CFj (T)^ on each volume. 
However, det ^[CFj (T)CFj (T -|-1)“^]^CFj (T)^ on each volume, which approaches det(C'^'l'C'®) with large T, has 
relatively large errors and fluctuations with no clear plateau and we fail to extract det((7^^C'®). This would be due to 
the smallness of the signals in smeared-point correlators. Meanwhile, det (^[Cfj^{T)Cfj^{T + 1)“^]^CFj^(T)^ shown 
in the middle panel in Fig. ^1 and det ([CFj^(T)]“^CFj (T)) shown in the right panel in Fig. which approach 
det(C'^'l'C'^) and det((C''^)“^C®) respectively, show relatively clear plateaus. Therefore we extract det(C'^'fC'^) 
and det((C'W)-iC'®) by the fits det(C'PtC'W) = det {T)CfJ^{T + l)-^f CfJ^(T)^ and det((C^)-1(7®) = 

det ([CFj^(T)]-icFj (T)) in the T range where they show plateaus and we finally obtain det((7^^(7®) as det(C'^'l'C®) = 
det((7^^(7'^) X det(((7^)“^(7®). In Fig. 1141 we show det(C^^C'®) obtained by the prescription shown above. The solid 
line denotes the best-fit curve by A\/V and the dashed line does the best-fit curve by AijV'^. The best fit parameters 
are A\ = 1.35 and Ai = 2.17, and x^/IVdf is 1.72 and 7.13 respectively. The volume dependence of det((7^'f(7®) seems 
not to be inconsistent with y. If we know the precise volume dependence of overlaps of wall operators, we may be 
discriminate the states using det((7^'l'C'^) or det(((7''^)“^(7®). 
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FIG. 14: Det(C'^^C'®) is plotted as the function of the lattice volume V. The solid line denotes the best-fit curve by AijV 
and the dashed line does the best-fit curve by . The best fit parameters are Ai = 1.35 and A 2 = 2.17, and /Ndt is 

1.72 and 7.13 respectively. These data behaves consistently in accordance with 1/17. 




